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First Measurement of r{D*^) 

CLEO Collaboration 
(February 5, 2001) 
Abstract 

We have made the first measurement of the D*^ width using 9/fb of e^e~ data col- 
lected near the T(4S) resonance by the CLEO II. V detector. Our method uses advanced 
tracking techniques and a reconstruction method that takes advantage of the small vertical 
size of the CESR beam spot to measure the energy release distribution from the D*^ — > 
D^TT^ decay. Our preliminary result is T{D*^) = 96 it 4 (Statistical) it 22 (Systematic) 
keV. 
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1 Introduction 



A measurement of r{D*~^) opens an important window on the non-perturbative strong physics 
involving heavy quarks. The basic framework of the theory is well understood, however, there 
is still much speculation - predictions for the width range from 15keV to 150 keV We know 
the D*^ width is dominated by strong decays since the measured electromagnetic transition 
rate is small, (1.68 ± 0.45)% 0]. The level splitting in the B sector is not large enough to 
allow real strong transitions. Therefore, a measurement of the width of the D*~^ gives unique 
information about the strong coupling constant in heavy-light systems. 

The total width of the D*~^ is the sum of the partial widths of the strong decays D*'^ D^n'^ 
and D*^ D^n^ and the radiative decay D*^ D'^'-f. We can write the width in terms of a 
single strong coupling, g, and an electromagnetic coupling, go+'y- 

r(D*+) = r(Z}°7r+) + r(D+7r°) + r(D+7) (i) 

2^' 3 ^ 9^ 3 ^ Q'^D+T 3 /ox 



using the isospin relationship 

9 = 9d+7t0 = goo-K+l v^, (3) 

where a is the fine structure constant, and the momenta are those for the indicated particle 
in the appropriate D*^ decay in its rest frame. The contribution of the electromagnetic decay 
can be neglected to note that the width of the D*^ only depends on (7 |]l|. 

Prior to this measurement, the D*^ width was limited to be less than 131 keV at the 90% 
confidence level by the ACCMOR collaboration |^. The limit was based on 110 signal events 
reconstructed in two decay channels with a background of 15% of the signal. This con- 
tribution describes a measurement of the D*^ width with the CLEO II. V detector where the 
signal, in excess of 13,000 events, is reconstructed through a single, well-measured sequence, 
D*'^ — i> 7r^Q„Z}°, — » K~7i'^. Consideration of charge conjugated modes are implied through- 
out this paper. The level of background under the signal is less than 3% in our loosest selection. 
Figure |1] compares the CLEO II. V signal with signal from the previously published limit. 

The challenge of measuring the width of the D*~^ is understanding the tracking system 
response function since the experimental resolution exceeds the width we are trying to measure. 
Candidates with mismeasured hits, errors in pattern recognition, and large angle Coloumb 
scattering are particularly dangerous because the signal shape they project is broad and the 
errors for these events can be underestimated, resulting in events that can easily influence 
the parameters of a Breit-Wigner fitting shape. We generically term such effects "tracking 
mishaps." A difficulty is that there is no physical calibration for this measurement. The 
ideal calibration mode would have a large cross-section, a width of zero, decay with a rather 
small energy release to three charged particles one of which has a much softer momentum 
distribution than the other two, and the other two would decay through a zero width resonance 
with a measurable flight distance. Such a mode would allow us to disentangle detector effects 
from the underlying width but no such mode exists. 
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D*-D Mass Difference for ACCMOR and CLEO II.V 
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Figure 1: The D*~^ - mass difference distribution obtained by ACCMOR (left)P] and in 
this worlc. 
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Therefore, to measure the width of the D*^ we depend on exhaustive comparisons between 
a GEANT based detector simulation and our data. We addressed the problem by selecting 
samples of candidate D*~^ decays using three strategies. 

First we produced the largest sample from data and simulation by imposing only basic 
tracking consistency requirements. We call this the nominal sample. 

Second we refine the nominal sample selecting candidates with the best measured tracks by 
making very tight cuts on tracking parameters. There is special emphasis on choosing those 
tracks that are well measured in our silicon vertex detector. This reduces our nominal sample 
by a factor of thirty and, according to our simulation, has negligible contribution from tracking 
mishaps. We call this the tracking selected sample. 

A third alternative is to select our data on specific kinematic properties of the D*~^ decay 
that minimize the dependence of the width of the D*^ on detector mismeasurements. The 
nominal sample size is reduced by a factor of three and a half and, again according to our 
simulation, the effect of tracking problems is reduced to negligible levels. We call this the 
kinematic selected sample. 

In all three samples the width is extracted with an unbinned maximum likelihood fit to the 
energy release distribution and compared with the simulation's generated value to determine a 
bias which is then applied to the data. These three different approaches yield consistent values 
for the width of the D*~^ giving us confidence that our simulation accurately models our data. 

2 CLEO Detector and Data Samples 

The CLEO detector has been described in detail elsewhere. All of the data used in this analysis 
are taken with the detector in its II. V configuration ||5|. This work mainly depends on the 
tracking system of the detector which consists of a three layer, double sided silicon strip detector, 
an intermediate ten layer drift chamber, and a large 51 layer helium-propane drift chamber. 
All three are in an axial magnetic field of 1.5 Tesla provided by a superconducting solenoid 
that contains the tracking region. The charged tracks are fit using a Kalman filter technique 
that takes into account energy loss as the tracks pass through the material of the beam pipe 
and detector p. 

The data were taken in symmetric e~^e~ collisions at a center of mass energy around 10 GeV 
with an integrated luminosity of 9.0/fb provided by the Cornell Electron-positron Storage Ring 
(CESR). The nominal sample follows the selection of D*^ n^^^D^ K^n'^ir^^^ candidates 
used in our — mixing analysis]^. 

Our reconstruction method takes advantage of the small CESR beam spot and the kine- 
matics and topology of the D*'^ — s> n^^^D^ — >• 7i^^^K~7r~^ decay chain. The K~ and 7r+ are 
required to form a common vertex. The resultant D° candidate momentum vector is then pro- 
jected back to the CESR luminous region to determine the production point. The CESR 
luminous region has a Gaussian width ~ 10 /im vertically and ~ 300 fim horizontally. It 
is well determined by an independent methodp|. This procedure determines an accurate 
production point for D°'s moving out of the horizontal plane; D^^s moving within 0.3 radi- 
ans of the horizontal plane are not considered. Then the vr^o^ track is refit constraining its 
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trajectory to intersect the production point. This improves the resolution on the energy 
release, Q = M{K~'!t~^7t^^^) — M{K~7r~^) — M(7r''"), by more than 30% over simply forming the 
appropriate invariant masses of the tracks. The improvement to resolution is essential to our 
measurement of the width of the D*^. Our resolution is shown in Figure |^ and is typically 150 
keV. The good agreement reflects that the kinematics and sources of errors on the tracks, such 
as number of hits and the effects of multiple scattering in detector material, agree between the 
data and the simulation. 

To further improve the quality of reconstruction in our sample, we apply some kinematic 
cuts to remove a small amount of misreconstructed background. Figure |^ shows the distribution 
of the momentum of the vr^o„ as a function of the D*~^ candidate momentum. We apply a cut 
at the kinematic boundary as shown in the figure. Figure ^ also shows the opening angle 6 
between the vr^ow ^^cl the candidate as a function of the D*"^ candidate momentum. We 
apply a cut of ^ < 38° /PD*[GeV] which is just beyond the kinematic limit to account for 
resolution smearing. We also require aq < 200 keV which removes the long tail in the error 
distribution. 

The tracking selected sample makes much more stringent cuts on the quality of the tracks 
used to identify the candidates. All tracks are required to have hits in both the r0 and z views 
in all three layers of the silicon strip detector as opposed to the nominal two silicon hits per 
view. None of these hits are allowed to be within 2 mm of a silicon wafer edge. The daughter 
tracks are required to have at least 38 of the possible 51 main drift chamber hits and seven of 
the ten intermediate drift chamber hits. The P^r degree of freedom of the fit to these two 
tracks are limited to less than 2 in each of the two drift chambers and 50 in the silicon strip 
detector. These selections are designed to remove tracks that have tracking mishaps or decay 
in flight. 

Since we have no calibration mode we compare the simulation and the data as a function 
of kinematic variables of the D*~^ decay. This will provide another test of the simulation's 
modeling of the data, and be the basis of our study of systematic uncertainties in the analysis. 
The most important kinematic variables are the "derivatives" which are defined by 



= mlo + m\ + 2{EdoE^+ - Pd^^P^+ cos 9) , (4) 



Pdo = Poo/Eno, (5) 



f)n E + 

^ - (/^Do - cos^), (7) 



dPoo M 
dQ _ Edo 



dP+ M 



(/?^+ -P no cos 9), (8) 



TT 



slow 



These derivatives test correlations among the basic kinematic variables, the and vr^o^ mo- 
menta and the opening angle, 9. We compare by dividing the Q distribution into ten slices 
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Gq , MeV 



Figure 2: Distribution of ctq, the uncertainty on Q as determined from propagating track 
fitting errors. The arrow indicates a selection discussed in the text. 
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Figure 3: The tt^^^ momentum on the left and the opening angle between _D° and n^^w 
right both versus the D*~^ momentum in the nominal data sample. 



in each of the kinematic variables and fitting the ten sub-distributions of Q to Gaussians. We 
display the width and mean of the ten fits as a function of each of the six kinematic variables 
in Figures ^ and |[ 

The quality of the width comparison (Figure ^ is excellent, with the simulation generated 
with an underlying r{D*~^) = 90 — 100 keV accurately predicting the data for all the kinematic 
variables. Even when generated with an underlying r{D*^) = keV the simulation accurately 
tracks the data's changes as the kinematic variables vary across their allowed range. 

The quality of the mean comparison (Figure ^) is not as good. The dependence of the 
mean of Q is not well modeled versus the vr^^^ momentum, dQ/dP^+ , and dQ/dPno by our 

slow 

simulation. We discuss the consequences of this imperfect modeling of the data in the section 
on systematic uncertainties below. 

Figure |^ shows the three derivatives plotted against each other in the data. Note that if 
we select dQ/dPoo and dQ/dP^+ both to be close to zero we minimize the dependence of 

slow 

Q on the basic kinematic variables Pno and P^+ , and thus minimize the contribution of the 

slow 

kinematic variables to the width of the Q distribution. With this selection we are more sensitive 
to the underlying width of the Q distribution rather than variations caused by any mismodeling 
of Q's dependence on the basic kinematics. The kinematic selection is defined by 



dQ 



dPoo 
dQ 



slow 



< 0.005, 

< 0.05. 
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Figure 4: Gaussian width of Q distribution versus kinematic parameters and derivatives. 
• - Data; o-Simulation with r^i.^ = 0; A-Simulation with Td*^ = 90 keV. 



10 



Data 



A -M.C. r=90keV 



^ 5.9 



A 

O) 

V 

c 

i 5.8 



a 



Jj I I I I I I I L 



_L 



J I I L_ 



5.9 



10 

, degree 



A 
V 

■2 
s 

75 



5.8 



d 



J I I I I I I I I I I I I L 



20 40 60 80 

dQ/dQ , MeV/rad 




K-slow 



, GeV 



0.1 
dQ/dP 



Tl-slow 



5.9 



A 

^5.85 



: t i 4 * * 



5.8 



J I I I 1_] I I I I I L_l_ 



5.9 



^5.85 



2.5 



3.5 



DO 



4.5 
GeV 



to 
s 



5.8 



t- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-0.01 -0.005 



0.005 0.01 
dQ/dP 



DO 



Figure 5: Gaussian mean of Q distribution versus kinematic parameters and derivatives. 
•-Data; A-Simulation with Td*+ = 90 keV. The horizontal hues show the average value 
of Q for the two samples. 
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Figure 6: Distributions of dQ/dPno versus dQ/dP^+ (top left), dQ/dO versus dQ/dP^+ 

slow slow 

(top right), and dQ/dO versus dQ/dP£,o (bottom) in the data. The dashed regions show the 
selection that defines the kinematic selected sample. 
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Table ^ summarizes the statistics in our three samples. 



3 Fit Description 

We assume that the intrinsic width of the D is nefflie: ible, T{D^) -C T{D*+), implying that the 
width of Q is simply a convolution of the shape given by the D*~^ width and the tracking system 
response function. Thus we consider the pairs of Q and ctq for D*~^ — >• ir^^^^D^ — >• K~7i~^7r^g^ 
where aq is given for each candidate by propagating the tracking errors in the kinematic fit of 
the charged tracks. We perform an unbinned maximum likelihood fit to the Q distribution. 

The shape of the observed Q distribution is assumed to be given by a signal that is a 
non-relativistic Breit-Wigner, 

s = r(g) 

{Q-Qo)' + {T{Q)/2r ^ ' 

with central value of Q, Qq. We also considered a relativistic Breit-Wigner as a model of 
the underlying signal shape, and found negligible changes in the fit parameters from the non- 
relativistic shape. The width of the signal Breit-Wigner is not a constant. It depends on Q 
and is given by 



where Tq is equivalent to r{D*'^), P and M are the measured candidate ir^^^ or momentum 
and D*^ mass in the D*~^ rest frame and Pq and Mq are the values computed using Qq. The 
effect of the mass term is negligible at our energy. 

For each candidate the signal shape is convolved with a resolution Gaussian with width 
ctq, determined by the tracking errors, as a model of our finite resolution. Figure ^ shows the 
distribution of aq for the data and the simulation. 

The fit also includes a background contribution with a fixed shape. The shape for the 
background is taken from fits to the background prediction of our simulation with a third order 
polynomial. The level of the background is allowed to float in our standard fit. The predicted 
background shape and fits are displayed in Figure 0. 

Figure |^ shows the Q distribution for our nominal data sample. Note that besides the well 
measured signal and the small, slowly varying background, there is also a small component 
centered on the signal with a large width. Therefore we allow a small fraction of the signal 
fmis to be parametrized by a single Gaussian resolution function of width amis- This shape is 
included in the fit to model the tracking mishaps which our simulation predicts to be at the 5% 
level in the nominal sample and negligible in both the tracking and kinematic selected samples. 
Typically we constrain the level of this contribution while allowing amis to float. 

We have many other parameters of the fit that can be varied or allowed to float for testing 
purposes. We can allow a scale factor on each candidate's ctq to model a systematic mistake 
in our tracking system caused, for example, by not properly accounting for the material of the 
detector. In our standard fits we only allow the normalization of the background to float, but 
we can either vary the shape as indicated by the simulation or allow the parameters of the 
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Evolution of the background shape 




Q,MeV 

Figure 7: Our simulation's prediction of the background for the three samples discussed in the 
text. Also shown are the fits to third order polynomials that are used in the fits to the data. 
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Figure 8: The Q distribution of the nominal sample in the data. 
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Table 1: Parameters of our fit to the Q distribution 



t*q rripf 


p c p T* 1 T-\ f 1 o n 


Tn 

u 


Breit Wigner width of D*+ Signal, r(D*+) 




Mean of D*+ Signal 


Ns 


Number of Signal Events 


fmis 


Fraction of Mismeasured Signal 




Resolution on Measured Q for Mismeasured Signal 




Number of Background Events 


k 


(Tq Scale factor, fixed to 1 




Coefficients of Background Polynomial, fixed from Simulation 



background polynomial to fioat as a measure of the small systematic uncertainty due to the 
background shape. 

Table |l] summarizes the parameters of our fit. Note that the ctq scale factor k and the 
background shape parameters -81,2,3 are fixed in our nominal fits. We minimize the likelihood 
function 

N 

L = 2(iV, + N,)-2J2 hg[NsS{Q,, ag,; Fq, Qo, frms, a^^s) + N,B{Qf, 51,2,3)], (14) 

i=l 

where 5" and B are respectively the signal and background shapes discussed above. 

The fitter has been extensively tested both numerically and with input from our full sim- 
ulation. We find that the fitter performs reliably giving normal distributions for the fioating 
parameters and their uncertainties. It also reproduces the input F (£>*+) from to 130 keV. Its 
behavior on each of the three data samples: nominal; tracking selected; and kinematic selected 
in the full simulation is discussed below. We note that if all the parameters are allowed to 
vary simultaneously there is strong correlation among the intrinsic width Fq, the fraction of 
mismeasured events fmis, and the ctq scale factor k, as one would expect. Thus our nominal fit 
holds k fixed, but in our systematic studies we either fix one of the three or provide a constraint 
with a contribution to the likelihood if the parameter varies from its nominal value. 

4 Fit Results 

As a preliminary test to fitting the data we run the complete analysis on a fully simulated 
sample that has about ten times the data statistics and is generated with a range of underlying 
r{D*~^) from to 130 keV. We do this for nominal, tracking, and kinematic selected samples. 
For the nominal sample we note that the fit is not stable if all the parameters are left to vary 
freely. We have found that if we constrain the fraction of mismeasured signal to 5.3 ± 0.5% 
as indicated by the simulation over the range of generated widths of the D*~^ then we get a 
stable result. This constraint makes the fit to the simulated nominal sample have no significant 
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Table 2: Results of the fits described in the text. The fit parameters are summarized in Table |I[ 
The uncertainties are statistical. 



Parameter 


Sample 


Nominal 


Tracking 


Kinematic 


To (keV) 


98.9 ±4.0 


106.0 ± 19.6 


108.1 ±5.9 


go (keV) 


5853 ± 2 


5854 ± 10 


5850 ± 4 




11207 ± 109 


353 ± 20 


3151 ±57 


fmis (%) 


5.3 ±0.5 


NA 


NA 


amis (keV) 


508 ± 39 


NA 


NA 


Nb 


289 ± 31 


15 ±7 


133 ±16 



offset between the generated and measured values for the width of the D*^. The tracking and 
kinematic selected samples have a negligible amount of mismeasured signal according to the 
simulation and in fits to these samples we fix f^is to zero. These simulated samples are also 
consistent with no offset between the generated and measured values for the width of the D*~^. 
We also note that in all three simulated samples there are no trends in the difference between 
measured and generated width as a function of the generated width; the offset is consistent 
with zero as a function of the generated width of the D*'^. Table |] summarizes this simulation 
study. We will apply these offsets to the fit value that we obtain from the data. 

Figures |lO|, and [ll| respectively display the fit to the nominal, tracking, and kinematic 
selected data samples. The results of the fits are summarized in Table |^. Correlations among 
the fioating parameters of the fit are negligible. Figure |T2| displays the likelihood as a function 
of the width of the D*^ for the fits to the three data samples. 

The agreement is excellent among the three fits, and when the offsets from Table |^ are 
applied we obtain 

Nominal Sample T{D*+) = 96.2 ± 4.0 keV, (15) 
Tracking Selected T{D*+) = 104 ± 20 keV, and (16) 
Kinematic Selected T{D*+) = 103.8 ± 5.9 keV. (17) 

The data sample and results are summarized in Table The uncertainties are only statistical. 
We discuss systematic uncertainties in the next section. 

5 Systematic Uncertainties 

We discuss the sources of systematic uncertainties on our measurements of the width of the 
D*+ in the order of their size. The most important contribution is the variation of the average 
Q as a function of the kinematic parameters of the D*~^ decay as shown in Figure ^ We 
estimate this uncertainty by repeating the fits described above in three bins for each of the 
six kinematic parameters and taking the uncertainty as the largest observed variation from the 
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Figure 9: Fit to nominal data sample. The different contributions to the fit are shown by 
different colors. 
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Figure 10: Fit to tracking selected data sample. The different contributions to the fit are 
shown by different colors. 
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Figure 11: Fit to kinematic selected data sample. The different contributions to the fit are 
shown by different colors. 
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Figure 12: Likelihood function versus measured D*^ width for the nominal (left), tracking 
(center), and kinematic (right) selected data samples. 



Table 3: Summary of our data sample, simulation biases, and fit results. 



Parameter 


Sample 


Nominal 


Tracking 


Kinematic 


Candidates 
Background Fraction (%) 

Tfit ~ Tgenerated (keV) 

Fit Fo (keV) 
D*+ Width (keV) 


11496 
2.51 ±0.27 
2.7±2.1 
98.9 ±4.0 
96.2 ±4.0 


368 
4.1 ± 1.9 
1.7 ±6.4 
106.0 ± 19.6 
104 ± 20 


3284 
4.05 ± 0.49 
4.3 ±3.1 
108.1 ±5.9 
103.8 ±5.9 
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nominal values in Table |^. We obtain uncertainties of ±16 and ±15 keV on r(Z)*+) and Qq 
respectively. 

The next most important contribution comes from any mismodeling of ctq's dependence on 
the kinematic parameters. We estimate this by varying our cut on aq from 75 to 400 from the 
nominal 200 keV and repeating our analysis with all parameters fixed except allowing the error 
scale factor k to vary freely. This indicates that the resolution is correct to ±4%, and we then 
repeat our standard analysis with k fixed at 0.96 and 1.04. We find uncertainties of ±11, ±9, 
and ±7 keV on r{D*~^) for the nominal, tracking, and kinematic selected sample. For Qo this 
uncertainty is negligible except in the tracking selected sample where it is ±4 keV. 

We take into account correlations among the less well measured parameters of the fit, such 
as k, fmisi and (Tmisi by fixing each parameter at ±lcr from their central fit values, repeating the 
fit, and adding in quadrature the variation in the width of the D*^ and Qq from their central 
values. We find uncertainties of ±8, ±9, and ±9 keV on the width of the D*^ for the nominal, 
tracking, and kinematic selected sample, and respectively ±3, ±4, and ±5 keV on Qq. 

We have studied in the simulation the sources of mismeasurement that give rise to the res- 
olution on the width of the D*^ by replacing the measured values with the generated values 
for various kinematic parameters of the decay products. We have then compared these uncer- 
tainties with analytic expressions for the uncertainties. The only source of resolution that we 
cannot account for in this way is a small distortion of the kinematics of the event caused by 
the algorithm used to reconstruct the origin point described above. This contributes an 
uncertainty ±4 keV on the width of the D*^ and ±2 keV on Qq. 

We consider uncertainties from the background shape by allowing the coefficients of the 
background polynomial to float. We observe changes on the width of ±4 keV for the nominal 
sample and ±2 keV for the tracking and kinematic selected samples. We have also released 
our kinematic selection cuts which causes the background to increase by a large factor. This 
causes a change which is small compared to allowing the coefficients of the background shape 
polynomial to float. Variations in the background have a negligible effect on Qq. 

Minor sources of uncertainty are from the width offsets derived from our simulation and 
given in Table |^, and our digitized data storage format which saves track parameters with a 
resolution of 1 keV and contributes an uncertainty of ±1 keV on the width of the D*^ and Qo- 

An extra and dominant source of uncertainty on Qq is the energy scale of our measurements. 
We are still evaluating the size of this contribution. 

Table ^ summarizes the systematic uncertainties on the width of the D*^ and Qq. 

6 Conclusion 

We have measured the width of the D*'^ by studying the distribution of the energy release in 
D*'^ D^'k'^ followed by K~tx^ decay. We have done this in three separate samples, one 

that is minimally selected, a second that reduces poorly measured tracks due to misassociated 
hits and non-Gaussian scatters in the detector material, and a third that takes advantage of 
the kinematics of the decay chain to reduce dependence on mismeasurements of kinematic 
parameters. The resolution on the energy release is well modeled by our simulation, with 
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Table 4: Systematic uncertainties on the width of the D*~^ and Qo 







Uncertainties in keV 










Sample 








Nominal 


Tracking 


Kinematic 


Source 


6T{D*+) 




6T{D*+) 


SQo 


6T{D*+) 


SQo 


Running of Q 


16 


15 


16 


15 


16 


15 


Mismodehng of ctq 


11 


< 1 


9 


4 


7 


< 1 


Fit Correlations 


8 


3 


9 


4 


9 


5 


Vertex Reconstruction 


4 


2 


4 


2 


4 


2 


Background Shape 


4 


< 1 


2 


< 1 


2 


< 1 


Offset Correction 


2 


NA 


6 


NA 


3 


NA 


Data Digitization 


1 


1 


1 


1 


1 


1 


Quadratic Sum 


22 


15 


22 


16 


20 


16 



agreement between the sources of the resolution as predicted by the simulation and analytic 
calculations. The largest sources of uncertainty are imperfect modeling of the dependence of 
the mean energy release on the kinematics of the decay chain, the simulation of the error 
on the energy release, and correlations among the parameters of the fit to the energy release 
distribution. With our estimate of the systematic uncertainties for each of the three samples 
being essentially the same we chose to report the result for the sample with the smallest 
statistical uncertainty, the minimally selected sample, and obtain 

r(D*+) = 96 ± 4 ± 22 keV, (18) 

where the first uncertainty is statistical and the second is systematic. We note that if we form 
an average value taking into account the statistical correlations among our three measures we 



get a result that is nearly identical with Equation |T8| since the average is dominated by the 
result with the smallest statistical uncertainty. This result is preliminary. 

This preliminary measurement is the first of the width of the D*~^, and our measurement 
corresponds to a strong couphng|]I|] 

^ = 17.9 ±0.3 ± 1.9. (19) 

This is consistent with theoretical predictions based on HQET and relativistic quark models, 
but higher than predictions based on QCD sum rules. 
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